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Background and Motivation 


r 


• Current LS-DYNA material models have been found to have 
limitations in the modeling of impact in composites 

- Existing models usually require significant a priori knowledge of 
damage and failure responses on the structural scale 

• Several years ago a new consortium was formed with the goal 
of creating a composite material model general enough to: 

- Model the wide range of material properties and architectures found 
in PMC’s 

- Recreate all of the behavior that can be found in material property 
tests (including tests that are not typically performed) 


Create a general material model for PMC composites which 
can predict impact results from mechanical properties, without 
relying on post impact correlation 
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• Developed model to be implemented into LS-DYNA as MAT_213. 

• Initial focus is on development of deformation model. 
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General Composite Model Requirements 


• Continuum Deformation/Damage Model with generalized, tabulated input, 
stress strain curve for non-damage related behavior (with limited or no 
curve fitting required by user) 

- Current models use point-wise properties that iead to curve fit approximations to 
actual material response 

- Tabulated input based on a well defined set of mechanical property tests leads to 
more accurate representations of actual material behavior 

• Input parameters based upon standard mechanical property tests - 
although alternate specimen test configurations or micro-mechanic 
analytical approaches producing virtual test results are acceptable 

• Effects of strain rate need to be accounted for in a flexible, unified 
manner accounting for anisotropy of rate effects. 

• Temperature dependency 

• Strain based damage and failure parameters 

• Shell and solid element implementations required (through thickness 
properties can be important) 

• Must be computationally extremely fast 
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General Approach 


r 


• Material nonlinearity in composites can be due to a 
combination of deformation and damage mechanisms 
acting independently or simultaneously 

- Most composite models assume that it is one or the other 

• The new model will allow both plasticity-like non-linearity 
and/or non-linearity caused by damage to be defined by 
the user 

- Tsai-Wu criteria (typically a failure surface) is used to define an 
orthotropic yield surface 

- Damage laws that model the orthotropic stiffness degradation will 
be defined by tabulated input 


Our goal is to create the capability to model general orthotropic behavior 
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Theoretical Formulation-Yield Surface 


Tsai-Wu failure criteria generalized to a yield function with 12 coefficients 
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Theoretical Formulation-Flow Surface 


Non-associative flow rule applied with 9 independent constants 
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Yield Surface Coefficient Evolution 


Current values of yield stresses required for evolution of yield surface 
and calculation of flow surface determined from 12 experimental (actual 
or virtual) stress-strain curves converted to be a function of effective 
plastic strain. 

A A f ^n(^") 





User provided load curves are 
true stress versus true plastic strain 





Internally stored is true stress 
versus effective plastic strain 
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Characterization of off-diagonai terms in yieid function 


Experimental and/or numerical variability can result in non-convex yield 
surface using default characterization method. 



Initial curves can lead to hyperbolic 
yield surface where elliptical surface is 
required 
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Off-diagonal terms adjusted to ensure 
convex yield surface and off-axis input curve 
adjusted accordingly 
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Characterization of Fiow Law Coefficients for Unidirectionai Composite 


• H11 , H12 and H13 identically equal to zero due to elastic response in 
composite longitudinal (1) direction. 

• H22 arbitrarily set equal to 1 by assuming in-plane transverse tension 
curve equal to effective stress-effective strain curve. 

• H33 equal to 1 due to transverse isotropy of unidirectional composite. 

• H23 equal to negative of plastic Poisson ratio v23. 

• H44 found by optimizing shear test or 45° off-axis test. 

• H66 equal to H44 due to transverse isotropy. 

• H55 found by optimizing “23” shear test or by using isotropic relation 
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Required Mechanical Property Tests to Characterize Composite Model 


High Strain Rate Tests 

Tension (1-direction) 
Tension (2-direction) 
Tension (3-direction) 
Compression (1-direction) 
Compression (2-direction) 
Compression (3-direction) 
Shear (1-2 plane) 


resulting test (data 


Damage & Failure Characterization Tests 

Tension (1 coupled w/2) 

Tension (2 coupled w/1) 

Tension (1 coupled w/3) 

Tension (3 coupled w/1) 

Tension (2 coupled w/3) 

Tension (3 coupled w/2) 
Compression (1 coupled w/2) 
Compression (2 coupled w/1) 
Compression (1 coupled w/3) 
Compression (3 coupled w/1) 
Compression (2 coupled w/3) 
Compression (3 coupled w/2) 


Material Characterization Tests 

ASTM Reference Test 

Tension (l-direction) 

ASTM D3039/D638 

Tension (2-direction) 

ASTM D3039/D638 

Tension (3-direction) 

ASTM D7291/D7291M-07 

Compression (l-direction) 

ASTM D3410/D3410M-03(2008) 

Compression (2-direction) 

ASTM D3410/D3410M-03(2008) 

Compression (3-direction) 

ASTM C365/D695,ASTM D7137 

Shear (1-2 plane) 

ASTM D5379 

Shear (2-3 plane) 

ASTM D5379 

Shear (1-3 plane) 

ASTM D5379 

Off-axis tension (45°, 1-2 plane) 

ASTM D3039/D638 

Off-axis tension (45°, 2-3 plane) 

ASTM C365/D695,ASTM D7137 

Off-axis tension (45°, 1-3 plane) 

ASTM D3039/D638 


• Not all materials, architectures, an(d 
(designs will require the full suite of tests 
for accurate pre(dictions (through 
thickness properties can be very 
important in impact response) 

• Mat 21 3 will fully accommo(date the 
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Overview of Numerical Implementation of Material Model 
Radial Return Based Approach 


• Step 1 : Compute elastic trial stress and yield surface coefficients 
from the current values of the yield stress. 

• Step 2 : Check if state is elastic by computing value of yield function. 

• Step 3 : If state is plastic (value of yield function greater than zero), 
start secant iteration to compute effective plastic strain and stress 
state that leads to yield function being equal to zero (convergence). 
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Verification of Materiai Modei 




In-Plane Shear 
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Stress (psi) Stress (psi) 


Laminate Level Verification of Material Model 


+/- 15® Tension (1-2 plane) 


250000 



0 0.005 0.01 0.015 0.02 

Strain 


+/- 45® Tension (1-2 plane) 

30000 



Strain 


+/- 30° Tension (1-2 plane) 



• Simulated curves show higher degree of 
nonlinearity compared to baseline most 
likely due to numerical method. 

• Discrepancy between experiment and 
analysis for [+!- 30-] curve most likely 
due to assumption that compression 
response equals tension response. 
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Proposed Revised Numerical Method to Account for Rotation of Yield 
Surface 


Yield Surface Evolution (1-2 plane) 
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Yield surface rotates with evolving plastic strain 
due to anisotropic yield function with yield stresses 
evolving in an anisotropic manner. 


Radial return algorithm adjusted to 
dynamically vary direction of plastic 
strain vector at each iteration and 
each time step. 
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Incorporating Strain Rate Sensitivity into Materiai Modei 



• Composite response can be extremely sensitive to strain rate. 

• Response in behavior at different rates of loading cannot be captured 
merely by scaling stress by rate. 

• Strain rate sensitivity will be incorporated into material model by use of 
tabulated input. 
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Stress (psi) 


Damage Coupling 

Triaxially Braided Composite Exampie 


• Damage created by strain in one direction can create damage (reduction 
in stiffness) in another, not correlated to the Poisson’s effect 

• Modeling of this effect will require an anisotropic coupled damage law 

- Most damage models are uncoupled 

• Testing performed at NASA GRC by Jon Salem and Nathan Wilmoth 


Stress V. Axial Strain: Tube 1035B 



0 - 0.0002 - 0.0004 - 0.0006 - 0.0008 - 0.001 - 0.0012 - 0.0014 

Axial Strain (in/in) 
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Preliminary Damage Modei 


Damage 
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Damage likely due to a 
combination of mechanisms. 

Most damage models assume only 
cracks perpendicular to load are 
present. 

Current model assumes 
“predamage” already exists due to 
transverse loads. 
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Required Mechanical Property Tests to Characterize Composite Model 


High Strain Rate Tests 

Tension (1-direction) 
Tension (2-direction) 
Tension (3-direction) 
Compression (1-direction) 
Compression (2-direction) 
Compression (3-direction) 
Shear (1-2 plane) 


resulting test (data 


Damage & Failure Characterization Tests 

Tension (1 coupled w/2) 

Tension (2 coupled w/1) 

Tension (1 coupled w/3) 

Tension (3 coupled w/1) 

Tension (2 coupled w/3) 

Tension (3 coupled w/2) 
Compression (1 coupled w/2) 
Compression (2 coupled w/1) 
Compression (1 coupled w/3) 
Compression (3 coupled w/1) 
Compression (2 coupled w/3) 
Compression (3 coupled w/2) 


Material Characterization Tests 

ASTM Reference Test 

Tension (l-direction) 

ASTM D3039/D638 

Tension (2-direction) 

ASTM D3039/D638 

Tension (3-direction) 

ASTM D7291/D7291M-07 

Compression (l-direction) 

ASTM D3410/D3410M-03(2008) 

Compression (2-direction) 

ASTM D3410/D3410M-03(2008) 

Compression (3-direction) 

ASTM C365/D695,ASTM D7137 

Shear (1-2 plane) 

ASTM D5379 

Shear (2-3 plane) 

ASTM D5379 

Shear (1-3 plane) 

ASTM D5379 

Off-axis tension (45°, 1-2 plane) 

ASTM D3039/D638 

Off-axis tension (45°, 2-3 plane) 

ASTM C365/D695,ASTM D7137 

Off-axis tension (45°, 1-3 plane) 

ASTM D3039/D638 


• Not all materials, architectures, an(d 
(designs will require the full suite of tests 
for accurate pre(dictions (through 
thickness properties can be very 
important in impact response) 

• Mat 21 3 will fully accommo(date the 
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Conclusions and Future Work 


r 


• New composite material modei MAT 213 being deveioped to provide 
improved predictive capabiiity for LS-DYNA simulations of composite 
impact. 

• Tsai-Wu composite failure model generalized to an orthotropic yield 
function. 

• Tabulated stress-strain curves used to track evolution of coefficients 
of yield function and stresses for flow law. 

• Characterization and numerical implementation of material model 
adjusted to account for issues related to orthotropic yield function. 

• Complementary damage model based on tabulated experimental 
input under development. 

• Strain rate and temperature effects being added to deformation 
model. 

• Further extensive sets of verification and validation studies planned. 
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